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Abstract 

<N 

We argue that conservation laws based on the local matter-only stress-energy-momentum tensor (characterized 
by energy and momentum per unit volume) cannot adequately explain a wide variety of even very simple physical 
phenomena because they fail to properly account for gravitational effects. We construct a general quasilocal conser- 
£ — . vation law based on the Brown and York total (matter plus gravity) stress-energy-momentum tensor (characterized 

by energy and momentum per unit area), and argue that it does properly account for gravitational effects. As a 
simple example of the explanatory power of this quasilocal approach, consider that, when we accelerate toward a 
freely-floating massive object, the kinetic energy of that object increases (relative to our frame). But how, exactly, 
does the object acquire this increasing kinetic energy? Using the energy form of our quasilocal conservation law, we 
can see precisely the actual mechanism by which the kinetic energy increases: It is due to a bona fide gravitational 
energy flux that is exactly analogous to the electromagnetic Poynting flux, and involves the general relativistic effect 
of frame dragging caused by the object's motion relative to us. 

> 

^ 1 Introduction and Summary 

m 

^ A conservation law ought to explain the change in some physical quantity contained inside a volume of space (e.g., 
^ total energy) in terms of related fluxes passing through the bounding surface of that volume. The standard approach 
(N to constructing conservation laws is based on the identity: V a (T ab ^ b ) = (V a T ab )^ b + T ab V^ b y where T ab is the 
• • matter stress-energy-momentum tensor, and ^ a is a vector that determines the type of conservation law, viz., energy, 
• »-j momentum, or angular momentum. In the context of general relativity, matter energy-momentum is locally covariantly 
^ conserved, i.e., V a T ab = 0, and the identity reduces to V a (T abl $ b ) = T ab V^ b y As is well known, the problem 
^ with this local conservation law is that the right hand side is, in general, not zero (or even a covariant divergence), 
resulting in a bulk term in the integrated conservation law that spoils what a conservation law ought to be. 

This problem is essentially gravitational in nature. There are several ways to see this: (1) The offending bulk term 
disappears when the spacetime has a suitable symmetry, i.e., admits a Killing vector, but dynamically interesting 
spacetimes (e.g., ones containing gravitational effects due to objects in motion) generically do not. The idea of relying 
on a spacetime symmetry to construct a conservation law is a throwback to pre-general relativity days, e.g., special 
relativity, where spacetime is maximally symmetric. The same goes for relying on asymptotic spacetime symmetries, 
where we are still in essentially a special relativistic mindset. Given that gravity is nontrivial spacetime geometry, 
an approach relying on spacetime symmetries cannot hope to properly incorporate gravitational effects in general. 
(2) The local conservation law above is homogeneous in T ab . In any matter-free region it is vacuous, even if that 
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region contains interesting gravitational physics, e.g., gravitational waves. The local conservation law is essentially 
blind to gravitational physics. We need a conservation law that is nontrivial even when T ab = 0. (3) Because of 
the equivalence principle, gravitational effects, e.g., gravitational energy, are not localizable, so we have no hope of 
capturing gravitational physics with a conservation law based on a local stress-energy-momentum tensor. For example, 
there is no such thing as a local gravitational energy density (energy per unit volume), that when integrated over a 
volume gives the total gravitational energy in that volume (see, e.g., §20.4 of reference [ 1 ]). (4) A bulk term in a local 
conservation law can be a symptom of the presence of fields that are not being accounted for in the stress-energy- 
momentum tensor. For example, in the standard Poynting theorem, the j ■ E bulk term is present because T ab excludes 
the charged matter field that is the source of the electromagnetic field, and represents an energy transfer mechanism 
between the electromagnetic field and the charged matter field. We contend that the bulk term in the local conservation 
law is, similarly, a result of T ab not properly accounting for the physics of the gravitational field. 

A solution to this problem is to move from local to quasilocal conservation laws, which can properly account 
for the gravitational physics. In this paper we construct a general quasilocal conservation law based not on the local 
matter stress-energy-momentum tensor, but on the Brown and York quasilocal total stress-energy-momentum tensor 
(matter plus gravity) 0. Here, "quasilocal" means that the differential conservation law is integrated not over the 
history of a volume of space, but over the history of the boundary of that volume. We focus on the case of energy 
conservation, and show that in the quasilocal approach, the quasilocal analogue of the offending T ab V f a ^b) bulk term 
becomes a surface flux term, which immediately solves the main problem mentioned in the opening paragraph above. 
Moreover, this surface flux term has two components: (I) The first component is a "stress times strain" term that 
can always be made to vanish by a suitable choice of frame, called a Rigid Quasilocal Frame (or RQF) 01 131 ■ The 
reader need not be familiar with RQFs to read this paper; it is sufficient to mention that, regarding point (1) in the 
previous paragraph, an RQF satisfies a certain "quasilocally projected" form of the timelike Killing vector condition 
for stationary spacetimes that allows us to move just far enough away from the spacetime symmetry mindset to include 
generic (i.e., non-stationary) spacetimes in conservation laws for energy, momentum, and angular momentum. (II) The 
second component is an "acceleration times momentum" term. This term is familiar from classical mechanics, and 
represents the rate at which the kinetic energy of an object increases due to one's acceleration toward it. Motivated by a 
simple equivalence principle argument, we show that this second term is actually a gravitational energy flux involving 
the general relativistic effect of frame dragging. We thus show precisely how quasilocal conservation laws resolve the 
bulk term problem in local conservation laws by properly accounting for the physics of the gravitational field. 

Our paper is organized as follows. In ^2]we introduce a very simple example of energy conservation in the context 
of special relativity, for the purpose of having a concrete example with which to illustrate the development of the 
general ideas. We consider a variant of Bell's spaceship paradox in which a box accelerates rigidly in a transverse, 
uniform electric field. Obviously, the electromagnetic energy inside the box increases, but how would co-moving 
observers explain this increase? Paradoxically, only half of the increasing energy comes from a net Poynting flux, and, 
according to the local energy conservation law, the other half comes from a bulk "acceleration times momentum" term 
integrated over the volume of the box. In £j3] we examine local conservation laws in general, with a particular focus on 
the role played by the T ab V \ a ^b) bulk term in our paradox example. This provides a point of comparison for in 
which we construct a general quasilocal conservation law and argue how it properly accounts for gravitational physics. 
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Figure 1: The trajectory of a rigidly accelerating cylinder 



We also apply the energy form of this quasilocal conservation law to the general relativistic version of our paradox 
example to concretely illustrate how, what we would normally think of as a bulk "acceleration times momentum" term, 
is actually a gravitational energy flux entering through the boundary of the box. In §[5] we present a complementary 
summary and conclusions, and argue that quasilocal conservation laws are necessary to understand more deeply a wide 
variety of phenomena, including the simple example of dropping an apple. 

2 An Apparent Paradox and its Resolution 
Paradox Outline 

Consider a right cylinder of length L and cross-sectional area A, whose axis is parallel to the Z-axis of an inertial 
reference frame in flat spacetime with Minkowski coordinates (cT, X, Y, Z). The cylinder is immersed in a constant, 
uniform electric field of magnitude E in the positive X-direction, and thus contains an electromagnetic energy density 
equal to E 2 /Sir. Now let us give the cylinder constant proper acceleration in the positive Z-direction so that the 
length (and volume) of the cylinder remain fixed for co-moving observers. It is well-known from Bell's spaceship 
paradox that, in order for the length of the cylinder to remain constant for co-moving observers, the top of the cylinder 
(represented by the hyperbola on the right in the diagram) must experience less proper acceleration, namely a' = 
a/(l + aL/c 2 ), compared to a, the proper acceleration experienced at the bottom of the cylinder (the hyperbola on the 
left) 0. The following two facts, also illustrated in the diagram, will be important to us: (1) The straight lines passing 
through the origin represent a natural choice for the surfaces of simultaneity for the co-moving observers: on any such 
simultaneity, all of the accelerating co-moving observers see the same instantaneous velocity, v, relative to observers 
at rest in the inertial frame. In other words, this is a "constant v" time foliation, with v monotonically increasing with 
time. (2) Co-moving observers at different positions along the length of the cylinder will see proper time flowing at 
different rates relative to one another. That is, between two co-moving simultaneities, if a proper time At elapses for 
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an observer at the bottom of the cylinder, a greater proper time, At' = (1 + aL/c 2 )Ar elapses for an observer at the 
top. 

With these facts in mind, let us consider the proper electric and magnetic fields seen by the co-moving observers. 
Let r (respectively, r') denote the proper time of observers at the bottom (respectively, top) of the cylinder, and (x, y, z) 
denote the natural choice of co-moving spatial Cartesian unit vectors. Since the observers are moving in a direction 
perpendicular to the electric field, they will see, in addition to a stronger electric field, a magnetic field; specifically, 
E = x*yE and B = —yf3 , yE, where /3 = v/c and 7 = — /3 2 . It is important to note that since these fields 

depend only on v, all observers on any given co-moving simultaneity will see the same instantaneous electric and 
magnetic fields. They will thus see the same proper Poynting vector, S = £^E X B = —zj^fi^ 2 E 2 , and the same 
proper electromagnetic energy density, u = ^(\E\ 2 + \B\ 2 ) = §^(1 + (3 2 )j 2 E 2 . Notice that, according to this 
expression for u, the total electromagnetic energy inside the cylinder is clearly increasing with time. We ask: How 
would the co-moving observers explain this increasing electromagnetic energy? One might expect that the change in 
energy between two co-moving simultaneities is just equal to the net Poynting flux into the cylinder over that time 
interval. As we will see, this is only half correct. 

Consider a pair of infinitesimally separated co-moving simultaneities labeled by proper time r and r + Ar as 
experienced by observers at the bottom of the cylinder. (A bit of thought shows that it does not matter whose proper 
time we use to parameterize the simultaneities.) Between these two simultaneities the total electromagnetic energy 
inside the cylinder changes by an amount 



where V = AL is the volume of the cylinder (which is constant for the co-moving observers). In this calculation 
we have made use of the relation dv/dr = a/7 2 for an observer at the bottom of the cylinder, who is experiencing 
constant proper acceleration a. 

Naively, we ought to be able to arrive at the same result by considering the net Poynting flux crossing the boundary 
of the cylinder. Before we calculate this, it is interesting to understand the mechanism by which the Poynting vector 
carries energy into the cylinder. As noted earlier, on any given co-moving simultaneity the relative velocity, v, is 
constant along the cylinder and, thus, the proper Poynting vector is the same at the top and bottom of the cylinder. 
Since no flux leaves or enters the sides of the cylinder, this suggests that the net Poynting flux is zero. Recall, however, 
that because of the acceleration, proper time advances more quickly at the top of the cylinder relative to the bottom. 
This results in a greater proper time-integrated flux entering the top of the cylinder compared to that leaving the bottom 
of the cylinder — in other words, a net accumulation of electromagnetic energy. Net electromagnetic energy enters the 
cylinder because of the time dilation effect associated with the acceleration. With the magnitude of the Poynting vector 
given above, we find 




(1) 




(2) 



Observe that the net Poynting flux, equation (|2]), accounts for only half of the change in electromagnetic energy inside 
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the cylinder, equation ([TJ. Apparently, we have a paradox. 



Paradox Resolution 

To understand the missing piece of this puzzle, let us first consider a simple problem in classical mechanics. Imagine 
an object with rest mass m from the point of view of a reference frame moving towards it with instantaneous speed 
v. In this reference frame, we would consider the object to have total energy E = \J m 2 c^ + c 2 p 2 , where p = ^ymv 
is the object's instantaneous momentum relative to the frame. Now, if our frame is accelerating at a proper rate a, it is 
easy to show that, in a proper time At, the object's total energy changes by an amount AE = apAr. Thus, the total 
energy of the object increases solely as a result of the acceleration of our frame relative to the object. 

The resolution to the paradox follows the same line of reasoning: as the cylinder accelerates relative to the existing 
momentum inside, there is an increase in energy proportional to this momentum times the acceleration. The momen- 
tum in this case belongs to the electromagnetic field contained in the cylinder, and is proportional to the Poynting 
vector: 

P=>> = ^ 2 (|V). (3) 

(As a consistency check, observe that as the electromagnetic energy, e = uV, inside the cylinder increases, so does 
the electromagnetic momentum, p, such that c 2 times the invariant mass, 

mc 2 = J e 2 - c 2 p 2 = — V. (4) 

remains constant, as expected, and equal to the electromagnetic energy inside the cylinder when v = 0.) 
Now, for the problem at hand, the contribution to the energy change due to this effect is just 

A£ Bu ik = apAr = 2/3 7 2 (~v) ~. (5) 

Notice that this energy change combined with the accumulation of Poynting flux due to the acceleration-induced time 
dilation along the length of the cylinder, Ai?p oynt i ng , now precisely matches that calculated from the energy density, 
i.e., AEpoynting + AE'suik = AE. Thus we have found the missing piece of the puzzle in the apparent paradox. It 
is not sufficient to consider the change in energy inside the cylinder based solely on the accumulation of Poynting 
flux across the surface; one must also take into account a "momentum times acceleration" term due to the frame 
accelerating relative to an existing momentum. The inclusion of this bulk term then resolves the apparent paradox 
but comes with a cost — it does away with the usual "change in energy equals net energy flux through the boundary " 
picture of a conservation law. 
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3 Local Conservation Laws 



Let us now look more generally at the notion of conservation laws in the context of both special and general relativity. 
Consider a smooth four-dimensional manifold, M., endowed with a Lorentzian spacetime metric, g ab , and associated 
covariant derivative operator, V a . In the presence of a non-zero matter stress-energy-momentum tensor, T ab , we can 
construct the identity 

V a (T ab V b ) = (V a T ab )* b + r o6 V (a * 6) . (6) 

This identity gives a differential conservation law for the current Q a = —T ab ^ b , whose physical interpretation depends 
on the choice of the weighting vector, ty a . Integrating both sides of this identity over a finite four-volume, V, we have 



- c J dZT ab u^ b = J dBT ab n a V b - j dV [(v a T ab ) V b + T ab (V (o tf 6) ) 



c 

£/-£; B V 



(7) 



On the left hand side, Ej and are the initial and final three-dimensional spatial volume "end caps" of V, with 
timelike future-directed unit normal vector \u^. On the right hand side, B is a three-dimensional timelike worldtube 
spanning the boundaries of the end caps, dT^ and <9£/, with spacelike outward-directed unit normal vector, n a , and 
induced Lorentzian three-metric ^y ab = g ab — n a n b . 

In general, the left hand side of equation ([7]) gives the change in a quantity contained inside a spatial volume (e.g., 
electromagnetic energy) over some time interval, while the first term on the right hand side corresponds to matter 
fluxes across the boundary of that spatial volume during that time interval (e.g., electromagnetic Poynting flux). This 
is the desired form of a conservation law. However, there are two additional terms on the right hand side that are 
four-dimensional bulk integrals. The first of these involves V a T ab , which is zero in general relativity; however in 
special relativity it will in general not vanish when T ab does not include all of the matter fields (e.g., when it does 
not include the electric four-current source of an electromagnetic field). The second term is present when \& a is not a 
Killing vector and, as it turns out, is crucial in resolving the apparent paradox above. To gain a better understanding 
of the terms in equation Q let us look at an example. 

Example: Electromagnetism 

Consider a general electromagnetic field, which can be decomposed as 

pab = 2 u [a E b\ + gafigc (8) 
C 

where u a is the four-velocity of a volume-filling, three-parameter family of observers who see proper electric and 
magnetic fields E a and B a , respectively, and e bc d = \u a e abc d is the spatial three-volume element orthogonal to the 
observers' worldlines. Denoting the spatial three-metric orthogonal to the observers' worldlines as h ab = gab-\-jsU a u b , 
the electromagnetic stress-energy-momentum tensor can be decomposed as 



T ab = JL [J_ u a u b < E 2 + B 2\ + 2^6) c ^ _ \gapb + fi « gb _ i^afe /g2 + ^2)1 \ 

4-7T [ 2c z c 2 ' J J 



(9) 
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Also, it follows from Maxwell's equations that 



V a T afe = j a F ab , (10) 

where j a is the electric four-current source of the electromagnetic field. 

For simplicity we will work with a congruence that has zero twist, i.e., one for which the four-velocity u a is 
hypersurface orthogonal, and assume that u a = uf, on and £j. This reduction in generality does not affect our 
conclusions — it just makes the formulas simpler and more transparent. Choosing as our weighting vector ^> a = \u a , 
equation (|7]) becomes an energy conservation equation, the individual terms of which are 

-T ab u^ b = u, (11) 
c 

T ab n a V b = - 1 n a S a , (12) 
c 



(V a T ab )^ b = j a E a , (13) 
T ab (V (a ¥ 6) ) = - c a a F a - T ab K ab . (14) 

Here u = ^ [E 2 + f? 2 ) and S a = ^e a bc E b B c are the proper electromagnetic energy density and Poynting vector, 
respectively, and j a is the electric four-current vector, as before. P a = ^S a is the proper electromagnetic momentum 
density, and T ab = —h a c h b d T cd is the spatial, three-dimensional Maxwell stress tensor. a a = u b V b u a is the observers' 
four-acceleration, and K ab = \h,°h b ^"^7 c Ud is the observers' spatial, three-dimensional strain rate tensor, measuring 



the expansion and shear of their congruence. Inserting equations (11 14 1 into equation (|7]) gives a general relativistic 



version of Poynting 's theorem for a hypersurface orthogonal congruence of observers: 

(15) 



d£u = -- / dBn a S a - I dV 



c 

Sy-Si B V 



j a E a + - o a P a - T ab K ab 
c 



Notice that the standard form of Poynting's theorem is recovered when u a is a timelike Killing vector, in which case the 



last two bulk terms on the right hand side, coming from equation ( 14 ), vanish. (The bulk term involving j a E a remains, 
however, and represents an energy transfer between the electromagnetic field and its electric four-current source.) 
When u a is not a Killing vector, however, these two bulk terms provide additional mechanisms for energy transfer. The 
a a F a term represents change in energy due to the frame accelerating relative to an existing electromagnetic momentum 
in the system (as was seen in resolving the apparent paradox of The T ab K ab term represents the Maxwell stress, 
T ab , doing work against the strain, K ab , of the three-parameter congruence. Since we would like to generalize beyond 
the case where u a is a Killing vector, both of these energy transfer mechanisms will be important. 
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Apparent Paradox Revisited 

We will now illustrate the use of the local conservation law, equation ([7]), to resolve the apparent paradox discussed in 
§[2] generalized slightly to allow for a time-dependent acceleration along the Z-axisQ Moreover, to facilitate compari- 
son with the general relativistic calculation in the next section (on quasilocal conservation laws), we will switch from 
a right circular cylinder of length L and cross-sectional area A to a round sphere of areal radius r. The reason for this 
switch is that the general relativistic calculation involves the extrinsic curvature of the boundary of our spatial volume, 
and the sharp corners of a cylinder introduce an unnecessary technical complication. 

Let us denote the Minkowski coordinates of an inertial reference frame in flat spacetime by X a = (X°, X 1 ) = 
(cT, X, Y, Z), I = 1,2,3. Let the embedding X a = £ a (r) define an accelerated, timelike worldline, Co, where r is 
proper time, and let e a (r) = \d£ a (r) / 'dr denote the unit vector tangent to Co- We will construct a three-parameter 
family of accelerated observers in the neighborhood of this worldline using the coordinate transformation 

X a (x)=e(r)+rr I (e, ( f>)e I a (T). (16) 

Here x a = (x , x l ,x % ) = (r, r, 9, 0) are coordinates adapted to the congruence of observers: the spherical coordinates 
(r, 9, <p) are parameters that label the observers' worldlines, and r is a time parameter along the worldlines which in 
general is not proper time, except for the observer at r = (i.e., there is a nontrivial lapse function, N, which will be 
given shortly). Also, r 1 (6, (j>) = (sin 9 cos 4>, sin 9 sin <fi, cos 9) are the usual direction cosines in a spherical coordinate 
system, and e 7 a (r) is a Fermi -Walker transported spatial triad defined on Co that is orthogonal to e a (r). 

Specializing to the case relevant to our apparent paradox, we let Co (i.e., the observer at r = 0) undergo proper 
acceleration a(r) along the X 3 (i.e., Z) axis, in which case our tetrad has the form: 

e a (r) = cosh a(r) 6% + sinh a(r) <5g, (17) 

e! a (r) = (18) 

e 2 a (r) = 5 a 2 , (19) 

e 3 a ( T ) = sinh a(r) 5% + cosh a(r) <$£, (20) 

where a(r) = - f r a(t) dt. A short calculation reveals that the metric in the observer-adapted coordinate system is 



9a/3 



( -c 2 N 2 

1 I - (21) 

\ r 2 Sy 



where E>ij = diagonal (1, sin 2 9) is the metric on the unitround sphere, and N(x) = l + ^a(r)r cos 9 is the lapse func- 
tion. The observers' four-velocity is defined as u a (x) = N~ l (x){d / dr) a = ce a (r), with resulting four-acceleration 
a a (x) = A r_1 (x)a(r)e 3 a (r). Comparing with the discussion of the paradox given in fj2j in the context of an accel- 
erating cylinder, the nontrivial lapse function, N(x), is the analogue of our earlier relation At' = (1 + aL/c 2 )AT, 



The analysis can be generalized to arbitrary acceleration; the limit to acceleration along the Z-axis is chosen purely for notational simplicity. 
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i.e., when a(r) > 0, proper time passes more quickly for observers above the plane 8 = tt/2 relative to those below 
it. Also, the magnitude of the observers' proper acceleration, A r_1 (x)a(r), is the analogue of our earlier relation 
a' = a/(l + aL/c 2 ), i.e., when a(r) > 0, observers above the plane 9 = -k/2 experience less proper accelera- 



tion compared to those below it. Finally, by inspection of the metric in equation (21 1, it is clear that the observers 
are accelerating rigidly, in the sense that the orthogonal distance between all nearest neighbor pairs of observers 
remains constant in time, as in the Bell spaceship paradox (despite the time-dependent acceleration). Thus, the three- 
dimensional strain rate tensor, K ab in equation ( [14] ), is zero. (This is a very special case. In general, it is not possible 
to find a volume-filling congruence of observers for which K ab = 0, a point we will return to in the next section.) 

As in £j2j we now introduce a constant, uniform electric field of magnitude E in the positive X-direction, 
perpendicular to the observers' motion: E a = E6f (and B a and j a vanish, at least in the region of the ob- 
servers' congruence). It follows from equation ([9]) that the electromagnetic stress-energy-momentum tensor is 
T ab = (E 2 /8tt) x diagonal(l, —1, 1, 1). Choosing \E ,a = \u a in equation 071) makes this an energy conservation 



equation. With = u a and n a = r 1 e 7 ° we find that the individual terms in equation (jT]) are 

-T^ b =(l + /3 2 ) 7 2 f^, (22) 

T ab n a ^ b = h 2 ^- cos 6, (23) 

4.7T 

V a T afe ) * b = 0, (24) 

T«» (V (a * 6) ) = -tf^jf, (25) 
where we have introduced the usual Lorentz transformation parameters: 7(7-) = cosh q(t) and /3(r) = tanh a(r) 



such that 7 = 1/y 1 — (3 2 . Observe that the right hand side of equation (25) is due entirely to the a a P a term in 



equation ( 14 1 since, as noted above, the collective rigidity of the observers' motion means that the strain rate tensor, 
Kab, vanishes in this case. In other words, this last contribution is entirely the result of a change in energy due to 
the frame accelerating with respect to the existing electromagnetic momentum in the system. Putting these results 
together, the conservation law reads 

f drr 2 d§ (l + /3 2 ) 7 2 — = f cdr r 2 dS/3 7 2 — cos9 [ 1 + \^(t)t cos^ + / cdr dr r 2 dS /3 7 2 — 
J 8-7T J 4-7T \ c l J J Air <r 

E / -S i B V 

(26) 

where the lapse function in the bulk term, T ab {y (a^b)) — see equation (25 1, has been canceled by that in the volume 
element, dV = cdr dr r 2 dS N. Note also that we have used dE = dr r 2 dS and dB = cdr r 2 dS N, where d§ = 
sin 9 d9 deft is the area element on the unit round sphere. 

In order to compare directly with the results in ^2j let the time interval be infinitesimal: tj — Ti = At. Using the 



relation Jp ((1 + /3 2 )7 2 ) = 4/37 2 |, the left hand side of equation (26 1 then integrates to precisely the same change 
in electromagnetic energy, AE, given in equation ([T]), except with the cylinder volume, V = AL, replaced with the 
sphere volume, V = Ivrr 3 . Also, the Ar in the equation now refers to the proper time elapsed for the observer at the 
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center of the sphere, instead of an observer at the bottom of the cylinder. The Poynting flux integral over B on the 



right hand side of equation (26 1 has two terms. The first term (proportional to cos 9) integrates to zero over the angles, 
which is analogous to the proper Poynting vector in S|2]being the same at the top and bottom of the cylinder such that, 
at lowest order, flux in equals flux out. At the next order in r, however, the nontrivial lapse function is responsible for 
a cos 2 6 term that does not integrate to zero. As in f]2j it is this non-isotropic time dilation that allows for a nonzero 
accumulation of Poynting flux. Evaluating this integral we find precisely the same change in electromagnetic energy, 
A^p oyn ti n g, given in equation Q, except with V and At reinterpreted as discussed above. Finally, the bulk integral 
over V on the right hand side of equation (26 1, which measures the change in energy due to the acceleration of the 
frame relative to the existing electromagnetic momentum in the system, evaluates to precisely the same change in 
electromagnetic energy, AEeuik, given in equation (|5]>, except with V and At again reinterpreted as discussed above. 

This example illustrates that any local conservation law constructed from the matter stress-energy-momentum 
tensor, T ab , as in equations (JHJ) and in general contains two terms responsible for the change in a physical quantity 
inside a volume of space: the first is a three-dimensional surface flux integral over B, as one might expect, and the 
second is a four-dimensional bulk integral over V. It is the addition of this bulk integral, which naively we did not 
expect, that resolves the apparent paradox introduced in ^2] But is such a bulk integral really necessary? In the 
next section we will see how to naturally convert this bulk integral into a surface flux integral, and at the same time 
generalize the conservation law in equations ([6]) and (|7]) to include gravitational effects. 



4 Quasilocal Conservation Laws 

We will now argue that the local conservation law given in equations ([6]) and (JTJ) is defective in two respects. First, 
even when matter is locally covariantly conserved, i.e., V a T ab = (which is always true in general relativity), the 
local conservation law contains a nontrivial bulk term when T ab is present and ^ a is not a Killing vector. Insofar as 
a generic spacetime does not admit any Killing vectors, when T ab is present this bulk term is generically present. But 
this bulk term violates what we usually think of as a conservation law, i.e., that the change in some physical quantity 
over a period of time is equal to some related surface flux integral over that period of time. Is there a natural and 
general way to convert this bulk term into a surface flux term? 

To help motivate this question, let us return for a moment to the "momentum times acceleration" bulk term that 
resolves the apparent paradox introduced in Jj2] We imagine being inside an accelerating box in flat spacetime that 
contains a freely-floating, massive object that appears to be accelerating toward us; the object's kinetic energy (relative 
to us) increases due to the acceleration of our frame. We ask: Where does the increasing kinetic energy come from? 
This might sound like a silly question — after all, energy is frame-dependent, and we are just changing the frame! 
However, the question is perhaps not so silly when we ask it in the context of the equivalence principle. Instead of 
being inside an accelerating box, we could imagine that the box is at rest in a uniform gravitational field, and that the 
object is experiencing an acceleration toward us due to the "force" of gravity. This "force" acting through a distance 
represents an energy transfer mechanism from the gravitational field energy to the kinetic energy of the object. So 
it might be possible to convert the "momentum times acceleration" bulk term into some kind of surface flux term 
representing gravitational energy entering the box from the outside. In the context of general relativity, i.e., when 
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we properly account for the frame dragging produced by the object in motion, we will see that this is exactly what 
happens. 

The second defect, which is related to the first, is that the local conservation law given in equations Q and ([7]) 
cannot properly account for gravitational effects since it is homogeneous in T ab . In any vacuum spacetime region 
where T ab = 0, the local conservation law has nothing to say. For example, we can imagine a vacuum region of 
space containing gravitational energy that is flowing in or out, to which the local conservation law is completely 
blind. This is, of course, to be expected, given that gravitational energy is not localizable. It seems that we need a 
term like V( a ^b)> which can be thought of as a measure of the presence of interesting gravitational physics, e.g., a 
non- stationary spacetime, and that this term should be coupled not to T ab , but rather some kind of quasilocal stress- 
energy-momentum tensor that represents both matter and gravity, so that it can be nonzero even when T ab = 0. We 
will presently construct a general quasilocal conservation law with precisely these properties, that will remove this 
second defect and, by its very construction, will also automatically remove the first defect, in an interesting and subtle 
way, exactly as anticipated in the previous paragraph. 

Let us consider an identity exactly analogous to equation ([6]), except defined in the three-dimensional spacetime 
of the timelike worldtube, B, i.e., the history of the boundary of a three-dimensional system. Then, for an arbitrary 
vector field, ip a , tangent to B, we have the identity 

D a (T§%) = (D a T^ b + n b D { M) , (27) 

where D a is the covariant derivative with respect to the three-metric, 7^ = g ab — n a n b , induced in B. In place of 
the four-dimensional matter stress-energy-momentum tensor, T ab , used in equation we have inserted the three- 
dimensional total (matter plus gravitational) stress-energy-momentum tensor, T^ b , defined by Brown and York 
This quasilocal stress-energy-momentum tensor is defined as Tg b = — - U ab , where U ab is the gravitational momen- 
tum canonically conjugate to the three-metric 7^ on B, and k = 8irG/c 4 . U ab , in turn, is equal to Q ab — @j ab , 
where @ ab = ^y a c \/ c n b is the exttinsic curvature of B. It will be useful to decompose the Brown- York quasilocal 
stress-energy-momentum tensor into energy, momentum, and stress surface densities, respectively, as: 



(28) 



£ = ^u a u b T® b [Energy/ Area], 

V a = -^cr b u c T^ c [Momentum/ Area], (29) 

S ab = a a c a b d T® [Force/Length], (30) 



where u a is the four-velocity of a two-parameter family of observers residing at the boundary of a spatial volume; 
i.e., the worldtube boundary, B, is the congruence of the integral curves of u a . The spatial two-metric, a ab = g ab — 
n a n b + \u a u b , projects tensors into the space orthogonal to both u a and n a , i.e., tangent to the spatial two-surface the 



observers reside on. Integrating equation ( 27 ) over a section of the worldtube, B, bounded by initial and final spacelike 
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slices, Si and Sf, we have 



dST$>u s a ^ b = I <U3 

Sf-Si B 



T ab n a ^ b 



(31) 



On the left hand side, dS is the surface area element on S{ and Sf, and ^-u^ is the timelike future-directed unit vector 
normal to Si and Sf (and tangent to B). 

It is interesting to compare this quasilocal conservation law to the local conservation law in equation (|7]). The left 
hand side of both equations gives the change in a physical quantity (e.g., energy) contained in a spatial three-volume. 



The first key difference is that the left hand side of equation ( 3 1 1 is quasilocal: it is an integral of a surface density 



(e.g., energy per unit area) over the two-surface boundary of the three-volume. Unlike the local volume density, the 
quasilocal surface density has no meaning by itself; only the quasilocal surface density integrated over a closed two- 
surface is physically meaningful. The second key difference is that the integrated surface density includes contributions 
from both matter and gravity (e.g., gravitational energy) 



The first term on the right hand side of equation ( 3 1 ) — the matter flux term — is identical to that in equa- 
tion (JTJ) (when \E' a = ip a is tangent to B), except its origin is very different. Unlike equation equation 
([27]) is purely geometrical, involving only the intrinsic and extrinsic geometry of B, with no reference to matter. 



The matter stress-energy-momentum tensor, T ab , enters equation ( 3 1 1 by first applying the Gauss-Codazzi identity 
(D a Tg b = — ^D a Il ab = —^n a G ab , where G ab is the Einstein tensor), and then applying the Einstein equation 
(G a b = KT a b). Thus, equation pTj ) forbids us from ignoring the fact that matter distorts the spacetime geometry; for 
example, it is meaningless to analyze an electromagnetic field in flat spacetime, as we did in £[3] Notice also that, 
because the boundary of a boundary is zero, there is actually no term in the quasilocal conservation law that is analo- 
gous to the matter flux term in equation ([7]). The matter flux term in equation ( |3Tj ) arises through an entirely different 
mechanism: it is actually a surface flux analogue of the j a E a bulk term in Poynting's theorem. The latter term arises 
when V a T ab ^ 0, which is analogous to D a T^ b ^ in the quasilocal conservation law. Specifically, V a T ab = j a F ab 
represents local energy-momentum transfer from the electromagnetic field to its sources, whereas D a Tg b = —n a T ab 
represents quasilocal energy-momentum transfer from the matter fields to the "system" contained in B. 



The last term on the right hand side of the quasilocal conservation law in equation ( 3 1 1 is analogous to the last 



term in equation yj. To make the comparison more detailed, let us look at energy conservation by taking 



Using the decomposition in equation (28 1 we find 

1 



pab 
l B 



D (aUb)) 



a a v a + s ab e ab , 



c 



(32) 



where a a = a a b a b is the projection of the observers' four-acceleration tangent to B, and 9 ab = \o~^a b ^ V c Ud is the 
strain rate tensor describing the time development (i.e., expansion and shear) of their congruence. The structure of 



equation (32i is identical to that of equation ( 14 1, except it is one dimension lower: surface, rather than volume, i.e., 
quasilocal, rather than local. By the very construction of the quasilocal conservation law, the volume densities a a ¥ a 



and —T ab K a h in equation (M_4b automatically become surface densities, a a V a and S ab 6 a b, respectively, i.e., surface 



abf 
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fluxes. This removes the first defect of local conservation laws discussed at the beginning of this section. 

Let us first focus on the quasilocal "momentum times acceleration" term, a a V a . Unlike the momentum volume 
density, P a , the momentum surface density, V a , is purely geometrical: it is a component of the extrinsic curvature of 
B closely related to a vector, cu a , that measures the precession rate of the observers' local "radial" vector, n a , relative 
to inertial gyroscopes. As shown in more detail in reference (U, the precise relation is: V a = e a b uj b , where 
tab = \tabcdU c n d is the two-dimensional spatial volume form orthogonal to both u a and n a . This allows us to write: 

= 7T7;£aba a Lu b , (33) 

07TG 

which is reminiscent of the normal component of the electromagnetic Poynting vector, h ■ x B. In fact, as we 
argue in more detail in reference [4] , we believe that a a V a is the exact, purely geometrical, operational expression for 
gravitational energy flux. So the "momentum times acceleration" bulk term, a a ¥ a , has really become a gravitational 
energy surface flux term, exactly as anticipated in the equivalence principle argument given near the beginning of 
this section. In the context of general relativity, the presence of momentum (matter or gravitational) flowing through 
the system causes the observers' local "radial" vector to precess relative to inertial gyroscopes (i.e., a frame dragging 
effect), and the vector cross product between this precession rate (the gravitational analogue of a magnetic field) and 
the observers' acceleration (the gravitational analogue of an electric field) corresponds to a flow of gravitational energy 
into the system. 

Of course this frame dragging mechanism of energy transfer would be difficult to confirm experimentally because 

2 

the frame dragging precession rate is smaller than the momentum density by the factor . Nevertheless, the existence 
of this mechanism is not surprising given that it bears a close resemblance to a similar mechanism in electrodynamics. 
In the context of special relativity, consider a charged particle accelerating in a uniform electric field (for simplicity 
we will ignore radiation reaction, which doesn't affect the point of our discussion). We ask: Where does the particle's 
increasing kinetic energy come from? The answer is obviously "the field." But what is the mechanism, exactly? 
Consider that the velocity field part of the magnetic field due to the particle's motion is localized near the particle, and 
circulates around its axis of motion. A moment's thought shows that the vector cross product between the external 
electric field (causing the particle's acceleration) and this magnetic velocity field (moving with the particle) points 
toward the axis of motion, and is maximum in the plane containing the particle. In other words, there is a Poynting 
vector representing a flow of energy from the field to the particle. Intuitively, one might guess that this answers 
our question. Surprisingly, it seems that a detailed analysis of this basic energy transfer mechanism in classical 
electrodynamics has been done only very recently [6]. The authors of this reference verify that our intuition is correct, 
at least in the ultra-relativistic limitE] 

To the extent that the linearized approximation to general relativity is very similar to electrodynamics, one might 
expect essentially the same energy transfer mechanism to occur in gravitational physics, when, e.g., we replace the 
massive, charged particle accelerating in a uniform electric field with a massive (charged or uncharged) particle accel- 

2 In the slow-motion limit, only 2/3 of the particle's kinetic energy is supplied by the Poynting flux; 1/3 is supplied by a peculiar "interfer- 
ence" effect between the external electric field and the velocity field part of the particle's electric field [6|. This subtlety is interesting, but it 
does not affect the spirit of our discussion. 
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erating in a uniform gravitational field (or, via the equivalence principle, a frame accelerating toward a freely-floating 
particle). In the context of the linearized, slow motion approximation to general relativity, sometimes called gravi- 
toelectromagnetism (for a review, see Q), precisely such a mechanism has been confirmed; see, e.g., (8j|9l. What's 
new in this paper is that this frame dragging energy transfer mechanism is confirmed in the full, nonlinear general 
theory of relativity with arbitrary matter and with no approximations; the gravitational energy flux, a a V a , is an exact, 
purely geometrical contribution to the energy flux that must be considered to properly explain a wide variety of phe- 
nomena, including phenomena as basic as a falling apple. We will compute an explicit example of this energy transfer 
mechanism in the next section. 

Let us now focus on the quasilocal "stress times strain" term, S ab 6 a }„ in equation (32). The boundary observers' 
strain rate tensor, 9 a b, is the two-dimensional analogue of the three-dimensional strain rate tensor, K a b, in equa- 



tion ( 14 1, for a volume-filling set of observers. In the latter case, it is obvious that such a three -parameter family of 
observers cannot, in general, move rigidly, i.e., move in such a way as to maintain constant radar-ranging distances 
between all nearest neighbor pairs of observers. This is because the condition K a j, = represents six differential 
constraints on three functions (the three independent components of the observers' four-velocity, u a ). In other words, 
the local conservation law generically contains a nontrivial bulk "stress times strain" term, — T ah K a b. 

The same is not true for a ?wo-parameter family of observers, because there the condition 6 a b = represents only 
three differential constraints on the same three functions. In reference 0] we argue that, in a generic spacetime, it 
is always possible to find a two-parameter congruence of integral curves of u a , comprising £>, that is expansion- and 
shear-free. Moreover, the degrees of freedom of motion left to such a congruence of boundary observers are precisely 
the same as those for rigid body motion in Newtonian space-time, viz., three translations and three rotations, each with 
arbitrary time dependence. In other words, in the quasilocal conservation law (but not in the local conservation law), 
it is always possible to choose a family of observers such that the "stress times strain" term vanishes. We have called 
such a family of observers a Rigid Quasilocal Frame, or RQF for short. We refer the reader to references |H for 
a more detailed discussion of the properties of RQFs. For our present purposes we need only appreciate that, since 
a congruence of RQF observers is expansion- and shear-free, there are no extraneous fluxes entering or leaving the 
system due merely to a change in the size or shape of the system boundary. So all of the interesting physics lies in 
the a a V a surface flux term discussed above. This is a key advantage to working with a quasilocal conservation law in 
general, and RQFs in particular. We will compute an explicit example of an RQF in the next section. 

In summary, the quasilocal conservation law introduced in this section removes two defects of the familiar local 
conservation law: (1) The undesirable bulk terms in the local conservation law automatically appear as surface flux 
terms in the quasilocal conservation law, and (2) unlike the local conservation law, the quasilocal conservation law 
accounts for gravitational effects. For example, the quasilocal energy conservation law contains an energy transfer 
mechanism related to frame dragging that is needed to understand in detail what is actually happening with regards to 
energy flow when, say, an apple falls. In the next section we will apply the quasilocal energy conservation law to a 
general relativistic analogue of the apparent paradox introduced in <32] 
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Quasilocal Approach to the Apparent Paradox 

To illustrate the explanatory power of the quasilocal approach let us now look at a general relativistic analogue of 
the apparent paradox introduced in £j2] We will work with a spacetime whose metric is g a b in coordinates X a = 
(X^^X 1 ) = (cT, X, Y, Z), I = 1,2,3, and which contains a covariantly constant (i.e., uniform) electric field of 
magnitude E in the positive A-direction iPlOl fTTTl : 

ds 2 = _ (1 + X 2 /L 2 ) c 2 dT 2 + + d * 2 2/L2) + (1 ^yl /L2) + (1 " Y 2 /L 2 ) dZ\ (34) 



where L = c 2 / yGE is the length scale of the geometry. It is important to note that the presence of the electric field is 
intimately connected to the spacetime curvature — we are no longer working in flat spacetime. This is readily verified 
by checking that there exist non-zero curvature invariants such as R a bcdR abcd = 8/L 4 . Alternatively, it is easy to see 
that this spacetime is a product of two surfaces of constant curvature ifTOlfTTI . 

Observe that g a f, reduces to the Minkowski metric in the X = 0, Y = plane. In this plane let us construct a 
timelike worldline, Co, parameterized by proper time, r, and defined by the embedding X a = £ a (r), that represents 
a fiducial observer undergoing proper acceleration a(r) along the Z-axis. As in along Co we define a tetrad, 
e a (r) and e 7 a (r), where e a (r) = \ d^ a (r)j dr denotes the unit vector tangent to Co, and the spatial triad, e 7 a (r), is 
Fermi-Walker transported along Cq. Because g a ^ equals the Minkowski metric on Co, we can use the same coordinate 



components for this tetrad as given in equations ( 17 1-( 20 1. Around Co we will construct a two-parameter family (two- 
sphere's worth) of accelerated observers, whose congruence of worldlines comprises a timelike worldtube, B. Our 
notation will follow that used in Q The spacetime metric, g a t, induces on B the spacelike outward-directed unit 
normal vector field, n a , the Lorentzian three-metric, j a b = g a b — n a n^, and the covariant derivative operator, D a . 
Also, we let u a denote the observers' four- velocity, which is tangent to B, and a a b = g a b — n a nb + ^u a Ub denote the 
observers' spatial two-mettic. 

In ^3] we constructed a three-parameter (volume-filling) family of accelerated observers in flat spacetime in the 
neighborhood of Co using the coordinate transformation given in equation ( [To] ). For arbitrary fixed r, this coordinate 
transformation actually defines an exact RQF, i.e., a two-parameter family (two-sphere's worth) of observers who, 
despite their time-dependent acceleration, maintain constant radar-ranging distance from their nearest neighbors — 



they are moving rigidly. Indeed, inspection of equation ( 2 1 1 shows that ni and Uj both vanish, so the observers reckon 
they are on a surface with constant spatial two-metric aij = gij = r 2 Sjj , which is a round sphere of areal radius r. 
Now because spacetime is locally flat, at lowest order in r we can start with the same coordinate transformation given 



in equation (16), but because of the spacetime curvature, we must add corrections at higher order in r to achieve the 
same RQF, i.e., to achieve = r 2 Sjj. Thus we begin with the ansatz: 

X a (x) = C(r) + rr'ie, 0)e/(r) + J [F 1 (t, 9, 0)e/(r) + F°(r, d, 0)e o a (r)] + O (£\ , (35) 

where, as before, x a = (x°,x ,x l ) = (r,r,9,(j)) are coordinates adapted to the observers: the angular coordinates 
x l = (9, <j>) are parameters that label the observers' worldlines; the radial coordinate, r, parameterizes the size of the 
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RQF; and the time coordinate, r, represents our choice of simultaneities on B. The three arbitrary functions F 7 (r, 9, <fi) 
allow us to perturb the observers' worldlines to satisfy the three conditions aij = r 2 Sij needed to achieve a round 
sphere RQF of areal radius r. (An RQF of generic shape and size need only satisfy the three weaker conditions, 
ddij/dr = 0.) The arbitrary function F°(t, 9, 0) allows us to perturb the choice of simultaneities on B to achieve 
other natural geometrical conditions, which will be discussed later. For the time being, we will keep the four functions 
F 1 and F° arbitrary. 



Choosing ip a = -u a in equation (31 1 makes this a quasilocal energy conservation equation. Using GRTensorll 



running under Maple, we compute the individual terms in equation ( 3 1 1 and find: 



- c dST% b u s a ^ b = r 2 dS + 7 2 ^ sm 2 0(5 cos 2 »|l)r|(C + $ f + * F )r|^xO (^j j , 

(36) 

dB T ab n a ^ h = cdrr 2 d§ ^ 2 — cos 9 j 1 + ^ a(r)r cos 9 + O ( ^ j j , (37) 



dBTg b (D [a ip b) ) =dB {^a a V a +S ab 9 a ^ = -cdr r 2 /3 7 2 ^ sin 2 9*^ {l + O (~) } + dBS ab 9 ah , 

(38) 

where dS = sin9d9d(j) is the area element on the unit round sphere, as before, and 7(7") and /3(t) are the same 



Lorentz transformation parameters we defined after equation (25 1. We have also inserted the relation L = c 2 /yGE 



where appropriate. Let us analyze each of these terms separately. 



Equation ( 36 1 represents the quasilocal energy density. The first term on the right hand side, at order 1/r inside 
the braces, is a vacuum energy. Insofar as it is constant in time, this term does not contribute to changes in the energy 
inside the sphere, and so can be ignored. The quantity C is a function of the angles (9, 4>), but not r. So for the same 
reason just cited, it can also be ignored. The quantity is a homogeneous function of the angular derivatives of 
the F 1 , which is a total derivative. It integrates to zero for any choice of the arbitrary functions F 1 , and so can also 
be ignored. The quantity is an inhomogeneous function of the angular derivatives of the F 1 , which is a linear 
combination of the three round sphere RQF conditions, aij — r 2 Sij = 0. It is zero when we choose a set of functions 
F 1 that satisfy these conditions (which we will do later), and so can also be ignored. The only nontrivial term is the 
second one. Insofar as only the integrated quasilocal energy density is physically meaningful, we can subtract off the 
spherical harmonics that integrate to zero, leaving the effective quasilocal energy density: 

>^i^r r * dS {^ r+ ^ <39) 

Equation ( [37] ) represents the quasilocal Poynting flux. The quantity in braces is (at least up to the order indicated) 
equal to the lapse function, N( x ) = 1 + ^a(r)r cos 9 + fjrl > which is contained in dB. Up to the lowest nontrivial 
order in r, this lapse function agrees with that found in f|3j as before, it is this non-isotropic time dilation that allows 
for a nonzero accumulation of Poynting flux in the sphere. Note that, up to the order indicated, the Poynting flux 



is independent of our choice of arbitrary functions F 1 and F°. Equation (38) represents the sum of the quasilocal 
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gravitational flux and the quasilocal "stress times strain" term — recall equation (32). When we impose the RQF 
rigidity conditions on the functions F 1 , the latter term vanishes (i.e., the strain rate tensor, 9 a b, vanishes). This leaves 
only the quasilocal gravitational flux term, a a V a , to consider. Note that, up to the order indicated, the gravitational 
flux is also independent of our choice of arbitrary functions F 1 and F°. 

Putting these results together, with the assumption that the RQF rigidity conditions can be satisfied (we'll demon- 



strate this below), the quasilocal energy conservation law in equation (31 ) reads, at order r , 



F 2 f ( F 2 / 1 \ F 2 a(Y)r 

r 2 dS 7 2 r= / cdTr 2 d§ \(3j 2 —cos9 1 + -^a(r)rcos# )+(3~f 2 — sin 2 9-^- )■ . (40) 

127T / 47T V c z / air c z 



This quasilocal conservation law is to be compared with the local conservation law in equation ( 26 1. As we did for the 
latter equation, we will take the time interval to be infinitesimal: r/ — 73 = At. Using the relation £pj 2 = 2[3-f 2 * 



[compare this with 4- ((1 + /3 2 )7 2 ) = 4/37 2 ^, which differs by a factor of 2], the left hand side of equation (40) 
then integrates to precisely the same change in electromagnetic energy, AE, given by the local conservation law, and 
also equation |TJ) (with the cylinder volume, V = AL, replaced with the sphere volume, V = |vrr 3 ). The quasilocal 



Poynting flux integral over B on the right hand side of equation (40 1 is, as mentioned earlier, identical to the local 



Poynting flux integral in equation (26 1, and so yields the same change in electromagnetic energy, A£p oyn ti ng , as 
found in that case, and also given in equation Q (with, again, the cylinder volume replaced with the sphere volume). 
We emphasize once more the role of the non-isotropic time dilation (the lapse function in parentheses) that allows for 
a nonzero accumulation of Poynting flux. 



The quasilocal gravitational flux integral over B on the far right hand side of equation (40 1 is, as discussed earlier, 
the surface flux version of the bulk "acceleration times momentum" energy transfer term in the local conservation 
law. The difference is that the local momentum volume density, P a , has been replaced with the quasilocal momentum 
surface density, V a , and the integration is not over a volume, but the boundary of that volume. We stress again that 
the actual mechanism of the energy transfer is a general relativistic effect, viz., the vector cross product between the 
boundary observers' acceleration and the precession rate of their gyroscopes due to the frame dragging caused by the 
electromagnetic momentum inside the system. All together, then, half of the energy flux is due to the electromagnetic 
field (with the flux entering mainly near the poles of the sphere, i.e., cos 2 9), and the other half is due to the gravita- 
tional field (with the flux entering mainly near the equator of the sphere, i.e., sin 2 9). Both contribute to a changing 
electromagnetic energy on the left hand side of the quasilocal conservation law. 

One might wonder how a flux of gravitational energy through the boundary of the system becomes electromagnetic 
energy inside the system. It seems that in general relativity there is no distinction between the two forms of energy. 



All forms of energy are equivalent. This is not surprising when we look at the metric in equation (34 1. The geometry is 
nontrivial and so presumably contains gravitational energy (in some nonlocal, or quasilocal sense), and the nontrivial 
geometry in turn is intimately connected to the electric field. They cannot be disentangled. 

All that remains in this section is to verify that we actually can solve the RQF rigidity conditions, viz., the three 
differential constraints o~ij = r 2 E>ij on the three functions F 1 (r, 9, (f>) which ensure that the observers' frame is a rigid 
round sphere of areal radius r. As argued in reference [4], such RQF solutions always exist in a generic spacetime; 
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moreover, they are unique up to motions of the RQF equivalent to those of rigid motion in Newtonian space-time, 
viz., six arbitrary time-dependent degrees of freedom: three translations and three rotations. Aiming for the simplest 
solution, we first set F 1 (r, 9, <p) and F 2 (t, 9, <j>) to zero, which eliminates translations of the RQF away from Co in 
the X and Y directions. We next demand that the RQF is not rotating, i.e., that the twist of the observers' congruence 
is zero, i.e., that the observers' four-velocity, u a , is hypersurface orthogonal (as a vector field in B). This means that, 
by adjusting F°(t, 9, (/)) appropriately, we can always find a time foliation of B (i.e., choose the observers' surfaces of 
simultaneity) such that u a is orthogonal to them, which is equivalent to the two conditions n, = 0. Using GRTensorll 
we can readily find a solution that satisfies all of these conditions, given by: 

F 3 (r, 9, 0) = 1 7 2 cos 9(3 - cos 2 9) - 1 cos 9(cos 2 9 + 3 sin 2 9 cos 2 <j>), (41) 
6 6 

F°{t,9,4>) = -\(3j 2 cos 9(3- cos 2 9). (42) 

Note that the <j) dependence in F 3 (t, 9,<p) is not surprising since the electric field, being parallel to the <p = plane, 
breaks the azimuthal symmetry — the electric field affects the geometry, and hence the coordinate embedding of the 
RQF. One can show that these results are in agreement with the more generally derived formulas found in reference 
0. 

5 Summary and Conclusions 

Using the standard, local way of constructing conservation laws — see equations (|6]> and (|7]) — we analyzed conservation 
of energy in the context of a simple example in special relativity, viz., a box, rigidly accelerating along the Z-axis, that 
is immersed in a transverse, uniform electric field. According to the local energy conservation law, the electromagnetic 
energy inside the box increases due to two separate mechanisms: (1) Half of the increasing energy is due to energy 
flowing in from outside the box via a Poynting flux. Interestingly, even though the instantaneous proper Poynting 
vector is uniform throughout the box (suggesting no net flux), there is a net proper time-integrated flux due to the 
acceleration-induced relative time dilation between observers at the top and bottom of the box. (2) The other half of 
the increasing energy is not due to energy flowing in from the outside; rather, it is a bulk effect due to the co-moving 
observers accelerating relative to the existing electromagnetic momentum inside the box, i.e., an "acceleration times 
momentum" energy transfer mechanism familiar from classical mechanics, integrated over the volume of the box. 

One might wonder if these two energy transfer mechanisms adequately explain what's happening. The answer 
is: No. First, in special relativity we assume that spacetime is flat, even though there is an electromagnetic field 
present. This precludes any possible general relativistic effects such as a flux of gravitational energy associated with 
the curvature of the spacetime caused by the electromagnetic field. Secondly, even in the context of general relativity, a 
local conservation law cannot properly capture all of the gravitational physics; for example, gravitational energy is not 
localizable — there is no such thing as a gravitational energy per unit volume, so a local conservation law cannot tell 
the whole story. Another way to see this problem is to notice that the standard local conservation law is homogeneous 
in the matter stress-energy-momentum tensor, and so has nothing to say in a matter-free region of space, even though 
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that region may contain dynamical curvature, e.g., gravitational waves. 

To address these shortcomings in general, and in particular see what's really happening with regards to the in- 
creasing energy inside the box, we constructed a quasilocal conservation law based on the Brown & York quasilocal 
total (matter plus gravitational) stress-energy-momentum tensor defined on the history of the boundary of a spatial 



volume — see equations (27) and (31 ). Using the energy form of this quasilocal conservation law, we analyzed the 
general relativistic analogue of the simple example described above. We found, again, that the electromagnetic energy 
inside the box increases due to two separate mechanisms. The first mechanism is identical in form to the Poynting 
flux mechanism described above [number (1)], but conceptually it has a completely different origin: it is actually 
analogous to a surface flux version of the j • E bulk term in the standard Poynting theorem, except instead of energy 
being transferred locally from the electromagnetic field to the four-current source of the field, it represents energy 
being transferred quasilocally from the matter fields to the "system" contained inside the box. The second mechanism 
is, at first sight, conceptually identical to the "acceleration times momentum" energy transfer mechanism described 
above [number (2)], except the volume integral has been converted to a surface flux integral over the boundary of 
the box. Going further, we argued that this surface flux is actually a gravitational energy flux exactly analogous to 
the electromagnetic Poynting flux, E x B, with E replaced by the acceleration of the co-moving observers on the 
boundary of the box, and B replaced by the precession rate of their gyroscopes due to the frame dragging caused by 
the electromagnetic momentum (in general, matter and gravitational momentum) flowing through the box. In both 
cases, now, energy is entering from outside the box via surface fluxes: half is an electromagnetic energy flux (entering 
the box through its top and bottom) and the other half is a gravitational energy flux (entering the box through its sides). 
Because an electromagnetic field is intimately intertwined with the spacetime curvature it produces, general relativity 
does not distinguish between electromagnetic and gravitational energy entering the box — both contribute on the same 
footing to the increasing electromagnetic energy inside the box. 

We can understand the second, gravitational energy flux mechanism intuitively as follows. Imagine being inside an 
accelerating box in empty space, which contains a freely-floating massive body that appears to be accelerating toward 
you. From your perspective, its kinetic energy is increasing. Where does the increasing kinetic energy come from? We 
could "explain" it as simply the "acceleration times momentum" energy transfer mechanism familiar from classical 
mechanics. Alternatively, we could invoke the equivalence principle and say that our box is at rest in a uniform 
gravitational field, and that the kinetic energy is coming from the gravitational field outside the box via some kind of 
gravitational energy flux passing through the boundary of the box. An analogous mechanism exists in the context of 
electrodynamics that explains, in detail, how an accelerating charge acquires kinetic energy from the external electric 
field causing the particle's acceleration [6]; moreover, the gravitational version of this mechanism has been shown 
to exist in the context of the linearized, weak field approximation to general relativity (8j|9l- In this paper we have 
established the existence of this very basic, but subtle mechanism in the full, nonlinear general theory of relativity, with 
arbitrary matter and no approximations. Our analysis made use of the concept of Rigid Quasilocal Frames (RQFs), 
discussed more fully in references El |4j, in order to properly isolate the relevant energy fluxes passing through the 
boundary of the box. 

In conclusion, we have explained the bulk "acceleration times momentum" energy transfer mechanism familiar 
in classical mechanics exactly in terms of a simple, operationally defined, purely geometrical, general relativistic 



19 



gravitational energy flux passing through the boundary of the volume in question. Naively, one might argue that since 
there is no G or c in "acceleration times momentum", this cannot be a general relativistic effect. But it is. It is based on 
frame dragging (the gravitational analogue of "B" in "E x B"), which is a general relativistic effect. We don't notice 
this mechanism in our day-to-day experiences because the typical gyroscopic precession rate vector due to a nearby 

2 

object in motion is very tiny; but it is precisely this vector, multiplied by the huge number that we identify as 
the "momentum" of the object [more precisely, the quasilocal momentum surface density, rotated by 90 degrees — see 



equation (33 1]. This general relativistic gravitational energy flux mechanism is what's really happening in the bulk 
"acceleration times momentum" energy transfer mechanism in classical mechanics. This deeper understanding would 
not be possible in the context of local conservation laws, and is a nice example of why we need quasiloc&l conservation 
laws. 
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